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ABSTRACT  ( Continue  on  reveree  aide  if  necessary  and  Identify  by  block  number) 

The  electrostatic  field  associated  with  the  rotation  of  Jupiter,  relative  to  the 
rest  frame  of  the  satellite  Io,  would  be  distorted  if  the  satellite  were  an 
electrical  conductor  (as  the  data  on  Jovian  decametric  radio  emission  may 
require).  An  idealized  two-dimensional  model  of  the  distorted  electric  field 
configuration,  in  the  limit  of  a perfectly  conducting  satellite  or  satellite 
ionosphere,  has  been  constructed.  This  model  has  been  used  to  trace  the 
adiabatic  guiding -center  trajectories  of  energetic  protons  and  electrons  — 
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0 ABSTRACT  (Continued) 

.^across  Jupiter's  magnetic -field  lines,  which  are  taken  as  rectilinear.  The 
'adiabatic  trajectories  of  very  low-energy  particles  (cold -plasma)  are  thus 
found  to  avoid  the  satellite  and  escape  absorption.  In  the  limit  of  very  high 
particle  energies  the  adiabatic  trajectories  are  undistorted,  and  absorption 
proceeds  as  if  lo  were  an  insulator.^,  The  interpolation  between  these  limits 
/Ts  monotonic  for  protons,  such  thaT  (for  example)  Io  sweeps  out  a drift  shell 
half  as  wide  as  the  satellite  for  first  invariants  M~  1 GeV/G  (energies 
W - 20  MeV).  The  situation  for  electrons  is  more  complicated,  and  we  find 
no  absorption  from  adiabatic  trajectories  at  M £ 46  GeV/G  (W  <30  MeV). 
Electrons  having  M > 50  GeV/G  are  typically  swept  out  of  drift  shells  wider 
than  the  satellite  itself  (by  a factor~  1 . 5,  for  example,  at  M ~ 400  GeV/G, 
'OV  90  MeV) . However,  electrons  can  impact  only  a portion  of  Io's  exposed 
liemTsphere  for  M~  50-200  GeV/G  (W~  30-60  MeV)T^Thus,  the  particle- 
absorbing characteristics  of  an  electrically  conducting  Jovian  satellite  are 
found  to  depend  on  both  the  species  and  the  energy  of  the  incident  particle, 
and  the  satellite's  particle  - absorbing  cross  section  differs  systematically 
from  its  geometric  cross  section. 
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CHARGED-PARTICLF  ABSORPTION  BY  IO 


The  five  inner  satellites  of  Jupiter  remain  within  the  magnetosphere 
at  all  times.  They  execute  nearly  circular  equatorial  orbits  about  the 
planet.  These  facts  have  led  Mead  and  Hess  (1973)  to  predict  a depletion 
of  ze  noma  gne  tic  ally  trapped  (charged-particle)  radiation  at  drift  shells 
that  intersect  the  satellites.  Subsequent  quantitative  analyses  of  particle 
absorption  by  Jovian  satellites  (Coroniti,  1974;  Hess,  Birmingham,  and 
Mead,  1974;  Thomsen  and  Goertz,  1975)  have  typically  treated  the 
satellites  as  electrical  insulators  that  do  not  otherwise  (apart  from 
absorption)  perturb  the  adiabatic  trajectory  of  a charged  particle.  Such 
a treatment  (which  may  well  be  valid  for  Europa,  Ganymede,  Callisto,  and 
Amalthea)  identifies  the  absorption  cross  section  with  the  geometrical 
cross  section. 

However,  the  well-known  correlation  between  the  orbital  phase  of 
the  satellite  lo  and  the  terrestrial  observation  of  Jovian  decametric  radio 
bursts  (Bigg,  1964)  has  led  Piddington  and  Drake  (1968)  to  propose  that  lo 
is  a good  electrical  conductor.  The  detailed  consequences  of  such  a model 
for  lo  have  been  investigated  by  Goldreich  and  Lynden-Bell  (1969).  Whether 
the  conductivity  results  mainly  from  Io's  interior  composition  (Goldreich 
and  Lynden-Bell,  1969),  from  its  surface  composition  (Matson,  Fanale, 
and  Johnson,  1974),  or  from  its  ionosphere  (Kliore,  Fjeldbo,  Seidel, 
Sweetnam,  Sesplaukis,  Woiceshyn,  and  Rasool,  1975)  is  uncertain,  but  is 
of  little  consequence  to  the  present  discussion. 
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The  point  is  that  a conducting  Io  must  inevitably  distort  (in  the 


vicinity  of  Io)  the  electrostatic  field  E^(  = - W)  associated  with  Jupiter's 
rotation,  and  must  thereby  alter  the  adiabatic  trajectories  of  zenomagneti cal  1 y 
trapped  particles  there.  The  purpose  of  the  present  work  is  to  estimate 
the  quantitative  effect  of  this  distortion  on  the  ability  of  Io  to  absorb  charged 
particles  from  the  Jovian  magnetosphere.  Our  intent  is  not  to  prejudge 
alternative  radio-burst  mechanisms  which  may  not  require  Io  to  be  a good 
conductor,  but  rather  to  derive  some  further  consequences  of  those 
mechanisms  that  do  require  it  (cf . Huba  and  Wu,  1976). 

In  order  to  make  the  present  analysis  tractable,  we  have  adopted 
an  idealized  rectangular  model  in  which 


= <»o  B0/rl 


and 


B(r ) = - z 


.Bo  (ao/ri 


)3[1  - 


(x/rj)]. 


(1) 


(2) 


where  aQ  (=  71372  km)  is  the  radius  of  Jupiter,  a^  ( = 1930  km)  is  the  radius 
of  Io  (including  a substantial  part  of  its  ionosphere),  and  r^  (=  421600  km)  is 
the  radius  of  Io1  s orbit  (Blanco  and  McCuskey,  1961).  Thus,  Io  and  its  flux 
tube  constitute  a cylinder  (of  radius  r = a^)  that  extends  to  x = ±oo  , and  the 
magnetic  field  B^is  everywhere  antiparallel  to  Jupiter's  rotation  axis 

A ^ 

(ft  = z).  The  coordinates  (x,  y,  z)  originate  at  the  center  of  Io.  It  proves 

MM  AMfc 

convenient  to  introduce  also  the  cylindrical  coordinates  ( r , cp)  such  that 
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x = r cos  and  y = r sin  ©.  The  parameter  (=  4.  192  G,  according  to 
the  spherieal-harmor.ic  analysis  of  Smitn,  Davis,  Jones,  Coleman, 

Colburn,  Dyai.  and  Sonett,  1975)  represents  the  equatorial  dipolar  mag- 
netic-field intensity  at  the  surface  of  Jupiter. 

The  geometry  corresponding  to  (1)  and  (2)  is  illustrated  in  Figure  la. 
lo  and  its  flux  tube  constitute  /by  assumption)  a perfectly  conducting  cylinder 
in  which  Jupiter’s  magnetic  field  is  "frozen  '.  The  surrounding  magneto- 
spheric  plasma  flows  (along  equipotentials  of  V)  around  the  flux  tube  as  if 
the  plasma  were  an  incompressible  fluid  streaming  without  viscosity  past 
a solid  cylinder  (Goldreich  and  Lynden-Bell,  1969).  Guiding -center  drift 
in  the  y direction  here  corresponds  to  azimuthal  drift  in  the  magnetosphere 
of  Jupiter.  Moreover,  the  absolute  contributions  of  electric  drift  and 
gradient  drift,  as  derived  from  (I)  and  (2),  correspond  at  the  orbit  of  I o 
to  those  obtaining  in  Jupiter's  magnetosphere.  Since  the  coordinates  in 

Figure  1 are  Iocentric,  the  angular  velocity  Q in  (1)  must  be  assigned  a 

-4  - 1 

magnitude  (=  1.  3467  x 10  sec  ) equal  to  the  apparent  rotation  rate  of 
Jupiter  in  the  rest  frame  of  Io  (Blanco  and  McCuskey,  1961;  Donivan  and 
Carr,  1969;  Mead  and  Hess,  1973). 

Since  the  conductivity  of  Io  and  its  flux  tube  is  regarded  as  infinite  in 
the  present  idealization,  there  is  no  perturbation  of  B in  lo's  wake  of  the  sort 

MM 

described  by  Goertz  and  Deift  (1973).  Moreover,  the  present  model 
envisions  no  intrinsic  magnetization  or  susceptibility  for  Io,  as  had  been 
contemplated  by  Burns  (1968).  The  flow  velocity  (Q  r^  ~ 60  km/sec)  of 
Jovian  plasma  relative  to  Io  is  smaller  than  the  Alfven  speed  for  a (hydrogen) 


I 


Figure  1.  Drift  paths  of  representative  particles  in  the  vicinity  of  Io: 
(a)  particles  having  M = 0;  (b)  protons,  M - 1 GeV/G; 

(c)  electrons.  M = 40  GeV/G;  (d)  electrons,  M = 80  GeV/G; 
(e)  electrons,  M - 200  GeV/G;  (f)  electrons,  M = 400  GeV/G 
Dashed  curves  represent  critical  (grazing)  trajectories  and 
separ atrices , Other  trajectories  are  identified  by  values  of 
x /a,  or  of  <p . . 

oo  i M 


nlasma  density  N < 7 x 1H'  cm"  ^ ; Io  should  have  no  standing  bow  shock 


m :.h>  s-  -,-e  . . rver.  1975  . The  flow  would  be  subsonic  only  for  a plasma 
temperature  T > 10  eV,  but  violation  of  this  latter  condition  could  only 
produce  some  electrostatic  noise.  It  would  not  seriously  alter  the 
character  of  our  simple  model.  We  also  ignore  the  possible  effects 
if  a r-hoto«lectr ; ~ sneath  iGurnett,  1972;  Hubbard,  Shawhan,  and  Joyce, 
T*'.v  which  might  be  present.  The  major  weakness  of  the  present  model 
is  our  two-dimensional  geometry,  which  makes  Io  s flux  tube  infinitely 
ior.g,  make?  Ic  s orbit  a straight  line,  and  suppresses  the  curvature 
drift  ir.  Jupiter  s magnetic  field.  However,  the  results  obtained  here 
should  as  reasonably  valid  for  charged  particles  that  mirror  near 
(ideally  at)  Jupiter's  magnetic  equator. 

For  a particle  with  a vanishing  component  of  velocity  along  B, 
conservation  of  the  first  adiabatic  invariant  M and  total  (kinetic  plus 
potential;  energy  W imply  \ drift  trajectory  given  (e . g.  , Schulz,  1972) 
bv  the  ecuation 


[l  - (2/m0c2)  MB(r2]1/2  + («q/m0c2)  v 


= [1  + (2/mQc2)  MBj]1^2, 


where  q is  the  particle :s  charge  and  B^  is  the  field  intensity  at  the  point 
(r,  'j)  - (aj.Cj)  where  the  trajectory  intersects  Io's  flux  tube.  Since  V = 0 
at  r = a(I  the  right-hand  side  of  (3)  represents  the  conserved  quantity 
I A (W/m_c  ).  The  equipotential  trajectories  shown  in  Figure  la  cor- 
respond  to  M - 0 and  are  species -invariant.  Trajectories  for  protons 


and  electrons  having  other  values  of  M are  shown  in  the  remaining 
panels  of  Figure  1. 


Proton  trajectories  (see  Figure  lb)  that  just  graze  lo's  flux  tube 
at  cp  = 0 and  cp  = n are  separated  by  a distance  La^  there.  These  critical 
trajectories  (dashed  curves)  attain  a separation  a x < 2a^  when  mapped 
to  y = ico,  There  are  no  stagnation  points  (intersecting  trajectories)  in 
the  flow  pattern  (near  lo,  at  least  ) for  protons  having  M >0  nor  for 
electrons  having  M>  205.  2 GeV/G.  Critical  trajectories  of  the  latter 
(dashed  curves.  Figure  If)  just  graze  lo's  flux  tube  at  cp  = 0 and  -p  = rr  but 
attain  a separation  Ax  > 2a  j when  mapped  to  y = x oc  . Individual 
trajectories  can  be  labeled  by  x_^  (the  value  of  x at  y = ±oo)  or  by  cp ^ (the 
value  of  c?  at  r = a^).  Our  interpretation  of  Figure  1 is  that  lo  absorbs 
zenomagne  tic  ally  trapped  particles  from  drift  shells  (L,  in  the  sense  of 
Stone,  1963)  having  a bandwidth  AL  = Ax^/a0,  and  tJiat  particles  on  drift 
shells  outside  this  band  avoid  absorption.  Representative  values  of 
Ax  /2a,  are  listed  in  Table  1.  The  symbol  W denotes  the  (kinetic)  energy 
of  a particle  on  the  trajectory  x = 0 for  given  M.  The  E x B drift  becomes 

Aft* 

negligible  in  the  limit  of  large  M,  and  so  the  ratio  Ax  /2a  approaches 

oc  1 

unity  for  all  charged-particle  species  in  this  limit. 


The  values  of  Ax^Q/2a1  listed  in  Table  1 refer  to  guiding -cente r 
trajectories,  and  each  might  be  augmented  by  the  quantity 


-/  - ,,  2 3 ...  2 _ 2 3 1/2 

p/aj  « (2  mQc  M/q  B^  aQ) 


(4) 
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in  order  to  account  for  the  absorption  of  particles  whose  guiding  centers 

pass  within  one  gyration  radius  ("p  ) of  Io  (r  = a^)  at  y = 0.  The  sum 

Ax,^  = Ax^  + 2 p"  is  monotonic  with  M for  protons,  but  not  for  electrons 

* 3 

(which  show  a relative  minimum  in  Ax  at  M ~2  x 10  GeV/G).  However 

JO 

2 5 

protons  having  M>  10  GeV/G  and  electrons  having  M > 3 x 10  GeV/G  are 

3 3 • 

found  to  execute  guiding -cente r drifts  (2rry  m^TjC/qB^aQ)  y>  0.  2a^ 

during  one  gyration  period.  Thus,  a particle  with  given  M can  (over  a 

certain  range  of  asymptotic  impact  parameters  x_^)  either  collide  with 

Io  or  avoid  Io,  depending  on  the  particle's  gyrophase.  In  other  words, 

the  absorption  process  is  somewhat  probabilistic  (rather  than  deterministic  ) 

for  certain  ranges  of  x^.  However,  the  adiabatic  approximation  invoked 

in  the  present  calculation  of  Ax^/2a^  should  be  valid  enough,  since  we  find 

that  f / a , <0.3  for  j (ax  /2a.  ) - 1 > 0.3  and  3p"  / r , « 1 for  1 (/  x / 2a , ) 

» l ~ cc  i ~ r i jo  l 

— 1 ] <sc  1 . The  distorted  electric  field  and  its  characteristic  scale  length 
a^  are  important  in  the  former  contingency,  but  not  in  the  latter. 

Stagnation  points  in  the  flow  pattern  for  electrons  having  M<  200 
GeV/G  generate  the  critical  trajectories  (dashed  curves)  in  Figures  lc-ld. 
For  M<  46.  38  GeV/G  the  stagnation  points  occur  only  at  x = 0,  with 

y = ± [1  - (3cM/ q y-  r \ a^=  ± Ay /2,  (5) 

and  generate  a separatrix  (dashed  curve.  Figure  lc)  that  precludes  the 

adiabatic  impact  of  such  radiation-belt  electrons  with  lo  (hence,  Ax  = 0 

oo 

in  Table  1 ).  The  symbol  ~ (h  [ 1 + (2  MBp/m^c^)  (a^/r^)^]^^)  denotes 
the  particle’s  ratio  of  relativistic  mass  to  rest  mass  at  x = 0.  The  symbol 
y denotes  the  same  ratio  in  general,  and  stagnation  points  occur  wherever 
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• 86  xl  0 0.  9983  16.0225  4.55x10^  1.0766  0.  3739 

.08xl04  0.9995  50.  6677  1.44xl03  1.0230  1.1825 

. 94  xl  O'*  0.9998  160.2254  4.56,.  103  1.0072  3.  7393 


q E = (M/\»VB.  Thus,  the  stagnation  points  approach  y = ±a,  (and  the 
separatrix  degenerates  to  the  circle  r = a^)  in  the  limit  M = 0 (as 
illustrated  in  Figure  la).  The  stagnation  points  specified  by  (5)  recede 
to  y - + x as  M *50.  70  GeV  / G,  a limit  in  which  the  present  rectangular  ized 
field  model  is  not  satisfactory. 

The  range  of  M values  (46.  38-59.  74  GeV / G)  spanning  the  com- 
plicated transition  from  the  topology  of  Figure  lc  to  that  of  Figure  Id 
coincides  roughly  with  the  range  over  which  the  present  model  is  too 
crude.  Since  the  gradient  and  electric  drifts  nearly  cancel  at  x/aj 
< 10  for  electrons  in  this  range,  a good  description  of  the  transitional 
drift  trajectories  would  require  a more  accurate  model  of  the  variation 
of  E and  B with  radial  distance  from  Jupiter.  We  have  formulated  such 

MM  MM 

a model  and  have  found  its  consequences  to  agree  well  (within  1%  for  all 
values  of  Ax^/2aj  given  in  Table  1,  although  comparison  was  not  feasible 
for  electrons  having  M = 50-200  GeV/G)  with  those  of  the  present  model 
specified  by  (1)  and  (2).  However,  the  more  realistic  model  is  too  com- 
plicated for  adequate  description  in  the  space  available  here,  and  it  does 
not  yield  simple  analytical  criteria  such  as  (5)  for  locating  stagnation 
points  in  the  flow  pattern.  Thus,  we  view  the  present  (simplified)  model 
as  the  more  instructive,  despite  its  failure  over  the  range  M » 46-60  GeV/G 
(corresponding  to  electron  energies  W 30-35  MeV  at  Io).  We  explore 
further  consequences  of  the  present  model  in  this  spirit. 

For  electrons  having  M > 60  GeV/G  the  stagnation  points  appear  at 
y = 0,  with 
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X 


- ± [(3cM/qYrj£.B,)  - 1 ]_1/Z  (6, 

where  y = [~^  - 3 (y^  - 1)  (x/r^)]^.  The  iteration  implicit  in  (6)  converges 
quite  rapidly,  and  the  stagnation  points  thus  specified  generate  separatrices 
of  the  type  illustrated  in  Figure  Id.  Electrons  of  this  class  can  impact 
only  a portion  of  lo's  surface,  i.  e.  , the  portion  outside  a pair  of  longitude 
intervals  facing  toward  and  away  from  Jupiter.  The  stagnation  points 
approach  x = ±aj  (exposing  lo's  entire  surface  to  electron  impact)  for 
M~  200  GeV/G  (as  illustrated  in  Figure  le).  Actually,  the  outer  stagnation 
point  reaches  x = + a^  for  M = 199.6  GeV/G,  and  the  inner  one  reaches 
x = - a^  for  M = 205.  2 GeV/G  (above  which  there  are  no  stagnation  points 
in  the  electron-drift  pattern  near  Io:  see  Figure  If). 

The  results  summarized  in  Table  1 correspond  to  the  various 
adiabatic  cases  described  above.  These  results  suggest  that  lo's 
particle-absorbing  cross  section  differs  systematically  from  its  geometric 
cross  section,  and  in  particular  that  Io  cannot  absorb  electrons  having 
W<  30  MeV  from  Jupiter's  radiation  belt.  This  last  conclusion  seems  to 
conflict  with  the  observational  data  of  Fillius  and  Mcllwain  (1974),  which 
show  depletions  of  the  160-keV  and  9-MeV  electron  fluxes  at  the  orbit  of  Io. 
One  might  hope  to  reconcile  this  conflict  by  appealing  to  collisional  diffusion 
of  electrons  as  they  execute  the  curved  portion  of  the  separatrix  (and  nearby 
trajectories)  in  Figure  lc.  However,  such  an  effect  turns  out  to  be  negligible; 
only  particles  on  adiabatic  trajectories  that  pass  within  1 km  of  the 
separatrix  could  conceivably  experience  a substantial  (>  10%)  probability 
of  diffusing  (by  collisions)  across  the  separatrix  and  (thus)  onto  trajectories 


I 
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Some  alternative  diffusion  mechanism  may  be  responsible 


for  transporting  such  electrons  (W  < 30  MeV)  from  a significant  swath  of 
drift  shells  across  the  separatrix  to  produce  the  observed  absorption 


profile,  but  it  is  difficult  to  be  quantitative  about  this. 


The  accepted  global  diffusion  mechanism  (Brice  and  McDonough,  1973) 

for  Jupiter's  magnetosphere  involves  fluctuating  electric  fields  generated 

by  fluctuations  (perpendicular  to  Bj  in  the  circulatory  motion  of  Jupiter's 

neutral  atmosphere  at  (Jovian)  ionospheric  altitudes.  Hess  et  al.  (1974) 

2 

estimate  the  resulting  radial-diffusion  coefficient  as  D ~ 300  km  /sec  at 

2 

the  orbit  of  Io,  and  Coroniti  (1974)  estimates  that  D~250  km  /sec  there. 
Spatial  diffusion  produces  a mean-square  displacement  ~ 2Dt  after  an  inter- 
action time  t.  The  time  required  for  an  electron  to  execute  the  exact 
separatrix  adiabatically  turns  out  to  be  infinite  because  the  particle  is 
stationary  at  the  points  where  qE  = (M/~)  VB,  as  given  by  (5).  However, 
a rough  estimate  for  nearby  paths  is  that 


(ria ^ / r j Q ) [ 1 - (3cM/q: 


s>r3/z 


This  represents  a fraction  Ay/4r^  of  the  particle's  full  drift  period 
(2n/Q^)  between  encounters  with  Io,  where  Ay/2  is  the  radius  of  the 

separatrix  at  x = 0,  as  given  by  (5).  Representative  values  of  Ay/2a^,  t, 

1/2  2 

and  (2Dt)  , as  derived  from  (5)  and  (7)  with  D-300  km  /sec,  are  given 

in  Table  2.  Thus,  for  example,  electrons  on  adiabatic  trajectories  that 

3 2 

pass  within  ~ 1 0 km  of  the  separatrix  at  M = 20  Ge V/G  (2Dt  0.  08  a j 

5 2 

3 x 10  km  ) should  have  a probability  > 10%  of  diffusing  across  the 
separatrix  and  being  absorbed  by  Io  on  a given  encounter.  Tabulations  of 
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the  full  drift  period  (Zn/Q^)  between  encounters  with  lo  and  the  cor- 
responding root-mean-square  displacement  (4n D/ft  ) ^ suggest  that 
such  an  electron,  if  not  absorbed  on  a given  encounter  with  lo,  has 
little  chance  of  being  absorbed  on  its  next  encounter  (see  Table  2). 

The  foregoing  considerations  (based  on  D~  300  km^/sec)  apply 
only  if  Jupiter's  radial-diffusion  coefficient  is  unmodified  in  the  vicinity 
of  lo.  This  seems  unlikely  in  view  of  the  effect  that  lo's  conductivity  has 
on  the  adiabatic  trajectories.  An  evaluation  of  lo's  local  effect  on  the 
magnitude  of  D would  be  interesting,  but  would  exceed  the  intended 
scope  of  the  present  work.  It  would  also  be  interesting  to  calculate  the 
electrostatic  potential  V(r)  surrounding  an  lo  with  only  finite  conductivity. 
Adiabatic  particle  trajectories  in  such  a model  would  presumably  lie 
between  those  calculated  above  and  the  simple  trajectories  postulated  by 
Mead  and  Hess  (1973)  for  the  case  of  a perfectly  insulating  lo.  However, 
the  present  work  (even  in  the  absence  of  such  refinements)  has  demon- 
strated that  the  conductivity  of  lo  entails  a major  reconsideration  of  the 
mechanism  by  which  the  satellite  absorbs  charged  particles  from  the 
radiation  belt  of  Jupiter.  Conversely,  if  lo  is  observed  to  absorb  charged 
particles  in  the  manner  of  an  insulating  obstacle,  then  one  must  reconsider 
the  electrodynamic  aspects  (Piddington  and  Drake,  1968;  Goldreich  and 
Lynden-Bell,  1969)  of  lo's  interaction  with  the  magnetosphere  of  Jupiter 
and  the  consequences  of  this  for  the  theory  of  Jovian  decametric  radio 
emission. 
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Table  2.  Dimensions  of  Separatrix,  Drift  and  Diffusion  Times,  and  Root-Mean-Square  Displacements 
for  Zenomagnetically  Trapped  Electrons  at  the  Orbit  of  Io  (D  = 300  km  /sec). 
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THE  IVAN  A.  GETTING  LABORATORIES 
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The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 

fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
phe  .ic  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
% and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 

sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
r propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 

lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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